A Prediction from the Type III See-saw Mechanism 
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A simple ansatz that is well-motivated by group-theoretical considerations is proposed in 
the context of the type III neutrino see-saw mechanism. It results in predictions for m s jvn^ 

7— i ; 

\ and m T /uib that relates these quantities to the masses and mixings of neutrinos. 

(N ' 

^ | Simple unified models based on SO (10) and related groups can lead to the so-called "type III 

\Q '. Tl 

see-saw mechanism" for neutrino masses jlj. In the most general case the type III mechanism leads 

to a light neutrino mass matrix given by the formula M v = -(M N H + H T M^){u/Q), where M N 

is the Dirac mass matrix of the neutrinos, H is a dimensionless complex three-by-three matrix and 

\ u/Q. is the ratio of a weak-scale vacuum expectation value to a GUT-scale vacuum expectation 

value (VEV). In a subsequent paper the type III see-saw mechanism was shown to have certain 

a: .„„ 



form of neutrino mass matrices 



In the simplest case, where a minimal set of Higgs fields breaks 



B — L, one has H = I and the type III see-saw formula takes the simple form 

M v = -(M N + Mj;)^. (1) 

The main problem in constructing predictive models of neutrino masses and mixings with the 
usual "type I" see-saw formula M v = —MnM^MJ^, is to relate the Majorana mass matrix 
of the right-handed neutrinos Mr, with its six complex parameters, to measurable quantities. 



There are very special mode 



there is such a relationship 



ip H 



s, such as the recently much studied "minimal 50(10) models", where 



(For an exhaustive list of references on "minimal 50(10) models" 



see 0.) And the study of leptogenesis may tell us something about the structure of Mr (although 
leptogenesis has only a single data point to work with). In general, however, the lack of information 
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about Mr is a problem for the predictivity of type I see-saw models. (The so-called "type II see-saw 
mechanism" 6] assumes the existence of S , f7(2)^-triplet Higgs fields with small VEVs that couple 
directly to vlvl- About the type II mechanism we have nothing to say in this paper.) 

What makes the simplest version of the type III formula, given in Eq. (1), so remarkable and 
appealing is that it does not involve the masses of the superheavy right-handed neutrinos at all. 
As a consequence, the simplest type III formula opens the possibility of constructing models of 
quark and lepton masses that are extremely predictive. In particular, in models based on 50(10) 
or other groups that unify an entire family within a single multiplet, the Dirac mass matrix of the 
neutrinos Mn is typically closely related by the grand-unification symmetries to the mass matrices 
(also of Dirac type, of course) of the up quarks, down quarks and charged leptons, which we will 
denote respectively as Mjj, Mp, and M^. It is therefore possible in many models (for examples, 
see [3, 0, Ei | ) to predict the matrix Mn from a knowledge of the masses and mixings of the quarks 
and the masses of the charged leptons. This would allow, if Eq. (1) holds, the complete prediction 
of the mass ratios and mixing angles of the neutrinos with no free parameters. 

In this paper we will not be so ambitious. We have not found so far a full three-family model 
that is as predictive as that and where all the predictions (or "postdictions" ) are consistent with 
experiment. Rather, as an illustration of the possibilities of the type III framework, we will 
present here a simple ansatz for the heavier two families that is well motivated by group-theoretical 
considerations. This ansatz leads to two interesting predictions that are consistent with present 
experimental data. Before presenting the ansatz, we very briefly review the type III see-saw 
mechanism and formula. 

In models based on 50(10), there are two ways that the right-handed neutrinos Nf {i = 



1,2,3) can get mass, either through a renormalizable term such as 16jl6jl26#, or through a 
higher-dimension effective operator such as 16jl6jl6#16# / 'Mqijt ■ The former allows automatic 
conservation of "matter parity" , whereas the latter makes do with smaller multiplets of Higgs fields. 
In the latter case, the effective d = 5 operator arises most simply from integrating out three or 
more 5O(10)-singlets, which we will denote by l a or 5 a , that have the couplings Fj a 16jl a 16# and 
(Ms)ablalb- If only the Standard-Model-singlet component of the 16# has a non-zero VEV, and 
we denote it by 0, ~ Mqut, then one has the familiar "double see-saw" mass matrix: 



^"neutrino — {yii ^Yj ; 5 Q ) 



( (Mtf)y W »■ \ 

{Ml)ij o F ib n 



AT A? 

^ F aj fl (M S )ab J 



N? 

\s b J 



(2) 
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By integrating out the superheavy fields N[ and S a , one obtains M v = —M^M R l Mjj, where 
M R = -(FQ)M- 1 (FVL) T . This is just the type I see-saw formula, with an effective Mr. 

Now, if we assume that the 5J7(2)^-doublet Higgs field contained in 16# also gets a non-zero 
VEV (and there is no fundamental reason why it should not), and we denote it by u, then the 
double see-saw mass matrix takes the form: 



^-"neutrino — iyii > 





T 

N)ij 



\ 



(M N )a F ib u ^ 

F ib n 

P (M S )ab J 



F T .{ 



( r, \ 



N? 

{ s b ) 



(3) 



In this case, it is easy to show that the effective mass matrix of the light neutrinos takes the form: 

M'„- -MxMT.'Mi - (.U v + M l v)-. (4) 



-MnM^mJ, - {M N + M^, 



where, as before, Mr = —(FQ)Mg l (FQ) T . The first term is the usual type I see-saw contribution, 
and the second term is the type III see-saw contribution. (The origin of the type III term can be 
simply understood as follows. One can eliminate the uS and 5^ entries in Eq. (3), i.e. the entries 
Fu and F T u, by doing a rotation of the (fj, Nf) basis by an angle 9 = tan# = u/O,. That reduces 
the matrix in Eq. (3) to the same form as Eq. (2), but with the zeros replaced by terms of the type 
III form.) Both the type I and the type III terms in Eq. (3) are formally of order M^/Mqut- 
However, since the elements of Mn are actually small compared to My/ because of small Yukawa 
couplings (except perhaps for the third family), and since M^ comes in quadratically in the type 
I term but only linearly in the type III term, one might expect the type III term to dominate 
for generic values of the parameters. Moreover, in the limit that the elements of Ms are small 
compared to the GUT scale, the type I contribution becomes small. As was pointed out in Q], 
that is a good limit for the purposes of enhancing leptogenesis. It is therefore plausible that one 
can neglect the type I term, and we shall do so. 

Now let us turn to the ansatz for the various Dirac mass matrices. Suppose that these have the 
form (neglecting the small masses of the first family) 



M v 




0a 
b 1 



mu, Mb 




0c 
d 1 



m D , 



(5) 



M N 




g 
h 1 



mu, M L 



^ ^ 
e 

1° / l ) 



m D , 



(6) 
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where the "texture zero" in the 22 elements can be enforced by an abelian family symmetry, either 
discrete or continuous. We will say more on this later. And further suppose that the entries satisfy 
the conditions 

a + b = g + h, c + d = e + f. (7) 

The relations given in Eq. (7) are not arbitrary, but follow from group-theory if the elements of 
the mass matrices come from no operators except of the following simple types: 

(1) WilBjlOn, 

(2) 16*16,120^, 

(3) WiWjlOH^H/MGUT. 

(4) leiiGjie'uieH/MGur, 

Eq. (7) is satisfied no matter how many operators there are of any of these types. Any operator 
of type (1) gives a = g, b = h, c = e, and d = /, thus satisfying Eq. (7). Any operator of type (2) 
gives contributions that are flavor-antisymmetric (since the 120 is in the antisymmetric product 
of two spinors). Consequently, it gives a + b = 0, c + d = 0, e + / = 0, and g + h = 0, thus also 
satisfying Eq. (7) in a trivial way. 

Any operator of type (3) gives contributions of the form fifjVf[aQ(f) + (3Q(f c )]. Here Q is 
that generator of 5O(10) to which the VEV of the adjoint Higgs field (45 h) is proportional; Q(f) 
is the value of this charge for the fermion / (= u, d, £~ , u) ; Vf = v u or v& depending on whether / is 
of the weak-isospin up or down type; and the coefficients a and (3 depend on the way the 5O(10) 
indices are contracted in the operator. Thus an operator of type (3) will give, for instance, c + d oc 
(a + /3)(Q(d) + Q(d c ))v d ande + / oc (a+(3)(Q(£~) + Q(£ + ))v d . Since the terms did$H d and £-£ + H d 
must be invariant under the charge Q, it follows that Q(d) + Q(d c ) = —Q(Hd) = Q(£~) + Q(£ + ), 
and so c + d = e + /, satisfying Eq. (7). In the same way it is easily seen that a + b = g + h. 

Finally, consider an operator of type (4). One of the spinor Higgs fields (say the unprimed one) 
gets a superlarge VEV that breaks 5O(10) down to SU(5). The effective operator that results is 
then of the form (ccl0j5j + /?5jl0j)5ij, where the coefficients depend on the contraction of 5O(10) 
indices in the original operator. This gives no contribution to a, b, g, and h, and gives contributions 
to the other paramaters of the form c = / and d = e (note the transposition between Mp and 
Ml). Again, such contributions satisfy Eq. (7). 

Simple low-dimension operators that could give contributions not satisfying Eq. (7) 
are 16jl6jl26# (if the SU(5) 45 contained in the 126# got a non-zero VEV), and 
1 6 j 1 6 j 16^16^ / Mqut ■ 
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One might ask why we do not include the effects of operators of even higher dimension, such 
as 16jl6jlOH45^/M^ [/T , which are not obviously smaller than the dimension-five operators that 
we included in our analysis, and which would not satisfy Eq. (7) in general. Such operators ought 
indeed to be present. However there are reasons that one might expect them to be small, as we 
now explain. Consider the operator I62I63IO h 45 h /Mgut, which will contribute to the 23 and 32 
elements in our illustrative model. Since these elements are somewhat small compared to the 33 el- 
ements, either the effective Yukawa couplings in this term are small or the ratio (45#) /Mgut must 
be, or both. This operator can arise from integrating out a pair of multiplets 16' + 16' that have 
GUT-scale mass, as follows. Suppose the terms al63l6'l0//, 6I62I6 45#, and M16'16'. Integrat- 
ing out the primed fields gives an effective operator abl6 2 16 3 10 H 45 H M- 1 [l + \b45 H /M) 2 ^ 1 / 2 . If 
b or (45#)/M are small, then the higher order operators are highly suppressed. This is not to say 
that operators of higher dimension must always be unimportant, but it is a plausible assumption 
easily implemented that they can be neglected. 

To return to the texture zero in the 22 elements, it could be enforced, for example, by a U(l) 
family symmetry under which the I63, 10h and 16' H are neutral; the I62 has charge +1; and the 
45#, 16h, and 120^/ have charge —1. 

Given the simplest type III form (Eq. (1)), and the ansatz of Eqs. (5), (6), and (7), one has 









(a + b) 
^0 (a + 6) 2 



\ 



umjj 



(8) 



\ 





M L = e m D , (9) 
\^0 (c + d-e) I j 

with the quark matrices Mjj and Mr, given by Eq. (5). Consequently, to the extent that we can 
ignore the first family, the five parameters a, b, c, d, and e determine the following mass ratios 
and mixings of the second and third families: m c /m t , m s /mb, V c b, in^/m T , 1712/1713 (the neutrino 
mass ratio), = smO atm , and m r /m;,. Therefore there are two predictions. (We assume all the 
parameters are real.) 

What we have done is use the values of the five quantities m c /m t , m^/m T , V c b, 1712/1713, and 
Qatm to solve for the five parameters a, b, c, d, and e. Then we have used the resulting values of 
those parameters to "predict" the values of m s /nib and m r /m{, at the GUT scale. For the first 
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three inputs (m c /mt, m^/rrir, and V c b), which are fairly well known, we have taken the central 
experimental values and run them up to the GUT scale, assuming low energy supersymmetry. The 
running depends significantly on the value of tan f3, and so we make a predictions for a particular 
set of values of tan (3 that span the interesting range: 2, 3, 10, 25, 40, and 57. The other two inputs 
(7772/7773 and Q a tm) come from neutrino oscillation experiments (see the reviews ^j, U|) and have 
rather large error bars. (For example, Q a tm = 45° db 6° at Mz-) We have assumed hierarchical 
spectrum for neutrino masses with mi < 0.007 eV. Under this assumption the RGE evolved values 
of 7712 / ^3 an d Oatm at the GUT scale remain within 3% _of their low-scale values even for large 



13] 



121,1131 •) 

Hence, we drop their running 



tan /?. (For relevant renormalization group equations see 
and allow these two inputs to vary within the experimentally allowed range and plot our predictions 
for m s /mb and 7771-/7775 as a function of them in Fig. ^ 

We take the experimental values of the quarks from Ref. , except for m s for which we use 



the results of lattice calculations as given in Ref. ^ | and double the error as suggested in Ref. 1£ ] . 
The values of the CKM angles and the charged lepton masses are taken from PDG 2004 17]. 
In presenting our results for m s /mb and m T /rrib in Fig. ^ we give the percentage by which the 
predicted GUT values differ from the RGE-evolved experimental central values. 

In doing the renormalization group running we assume that all the sparticles have mass of 1 
TeV. From Mz to 1 TeV, the running is done at one loop, assuming the Standard Model with 
two Higgs doublets. From 1 TeV to the GUT scale (taken to be 2 x 10 16 GeV) we do a two-loop 
running assuming the MSSM. The gauge coupling constants are taken from PDG 2004 17]. We 
present one example of the RGE evolution in Table HJ 

It should be noted that, even with the assumption that we are making that the parameters a, b, c, 
d, and e are real, there are discrete ambiguities of the relative signs of these parameters. (The overall 
sign does not matter.) The choice that gives by far the best fits is (a, 6, c, d, e) = db(— , +, +, — , — ). 
A typical set of values is a ~ -0.00455, b ~ +0.9, c ~ +0.04, d ~ -0.45, e ~ -0.55, and 
/ = c + d-e~+0.14. 

Note that the value of a is very small. It is this that accounts for the smallness of m c /rrit. One 
way that a might be small naturally (i.e. without fine-tuning) using only the set of operators that 
satisfy Eq. (7) is by means of an operator of the form I62I63IO/745// /Mqut, where (45//) oc Q = 
ha + e(.B — L), e < 1, where I3R and B — L are the familiar S'O(IO) generators the diagonal 
generator of SU(2)ji in the Pati-Salam subgroup, and B — L the baryon minus lepton number), 
and where the fields are contracted in such a way that this generator Q acts on the field 162- (This 
would happen, for instance if the effective operator came from integrating out a 16' + 16' having 
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FIG. 1: Plots of constant values of m s /mb (solid lines) and m T /mh (dots) in the tan(9 atm -m2/m3 plane 
for tan (3 = 2, 3, 10, 30, 45, and 57. The percents represent the departure of the fitted ratio from the central 
value obtained from the RGE evolution. In the case of the m T /nib ratio the uncertainty in the b and r mass 
allows this departure to be within ± 8% range. 

the couplings 16216 45//, 163l6'10//, and M16'16'.) This operator would give off-diagonal mass 
terms for the up quarks proportional to Q(u)u^U2 + Q(u c )u2U^. Since, I^r(u) = 0, this would give 
a/b = 0(e). 

The values of m s /m b that we predict are satisfyingly close to the experimental (lattice) results. 
A couple of things should be noted in this regard. First, it was long thought that the Georgi- 
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TABLE I: Input values at the M z scale vs. the GUT scale values {M GUT = 2x 10 16 GeV) for tan/3 = 45. 
In the fermion case we use indicated errors at the Mz scale to extract corresponding errors for individual 
fermions at the GUT scale. No correlation is taken into account. 





ii = M z 


[l — Mqut 


tan/3(/i) 


45.00 


35.36 


Vu(ji) (GeV) 


174.05 


117.33 


v d (p) (GeV) 


3.87 


3.32 


m u (/i) (MeV) 


9 oq+0.42 
z - oo -0.45 




m c (fi) (MeV) 


677+1?; 


212+H; 


mt(ji) (GeV) 


181+11 


93+??; 


m d (p) (MeV) 


4.69l°;«° 


i r 9 +0.19 


m.{n) (MeV) 


53.8™ 


17 4++ 3 


m 6 (M) (GeV) 


3-00t°il 


1 <U+0-08 

J-^-o.o? 


m e (n) (MeV) 


n ASRSA797+ - 00000014 


35fi20421+ 000000010 

U.OODZU^Z1_q .00000010 


mpQj,) (MeV) 


102 751 4n+ 00033 


7 5 20781 + 00024 

'°- zu ' ol -0.00024 


mr(/i) (GeV) 


1 74669+ 00030 

l.(*OO»_ ,00027 


1 /I ci 1 1 +0.00032 
1.1Olll_ . 00029 


|(VcKAf)l2(M)| 


0.2200 


0.2199 


\(V C km)i3(h)\ 


0.00367 


0.00300 


\{V C km)2?,{p)\ 


0.0413 


0.0337 



Jarlskog jjg] relation (rris/m^GUT = \{m^/rnr)GUT gave a good fit to the data in SUSY GUT 
models. However, the recent lattice calculations have given results for m s that are typically only 
about 0.6 times the typical values that had been obtained by previous methods. Because of that, 
many models which were constructed in the past to give the Georgi-Jarlskog result, would be off 
from the current central experimental/lattice results for m s /m{, by about +60%. That compares 
to the values we are getting, which agree with the current central value of m s /mb for some of the 
allowed (7712/1713) — (0 a tm) parameter space, and are within 20% for a large part of that space. 

A second point is that inclusion of the first family is likely to push up the predicted value of 
m s /mb by about 5%. The reason is that empirically the relation for the Cabbibo angle 9c — 

n 

^/md/m s is known to work very well 19]. As is well-known, this formula arises naturally if the 11 
element of the down quark mass matrix vanishes and the 12 and 21 elements are approximately 
equal l^jj. But then diagonalizing the 12 block of the down quark mass matrix will push up the 
value of the 22 element by a factor of (1 + |m^/m s |). 

In any event, we see that further improvement in the measurement of the 6 a tm, Sm^ tm , 5m 2 sol , 
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and the lattice results for m s , together with an eventual determination of tan/3 will allow our 
simple ansatz, given in Eq. (7), to be tested. 
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